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Abstract 

In this paper we critically discuss the so far performed attempts 
aimed at the detection of the general relativistic gravitomagnetic Lense- 
Thirring effect in the gravitational field of the Earth with the existing 
LAGEOS satellites. In the latest reported measurement of the grav- 
itomagnetic shift with the nodes of the LAGEOS satellites and the 
2nd generation GRACE-only EIGEN-GRACE02S Earth gravity model 
over an observational time span of 11 years a 5-10% total accuracy is 
claimed at l-3cr, respectively. We will show that, instead, it might 
be 15-45% (1-3<t) if the impact of the secular variations of the even 
zonal harmonics is considered. Possible strategies in order both to 
make more robust and reliable the tests with the node-only LAGEOS- 
LAGEOS II combination used and to overcome the problems affecting 
it with other alternative combinations are presented. 
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1 Introduction 



Recent years have seen increasing efforts aimed to directly detecting various 
phenomena connected to the general relativistic gravitomagnetic field jU 
IH 03 of the rotating Earth. 

The extraordinarily sophisticated and expensive Gravity Probe B (GP- 
B) mission 17] has been launched in April 2004; it is aimed at the detection 
of the gravitomagnetic precession of the spins [S] of four superconducting 
gyroscopes carried onboard at a claimed accuracy of 1% or better. 

The Lense-Thirring effect on the orbital motion of a test particle jH] 
could be measured by analyzing the orbital data of certain Earth artificial 
satellites with the Satellite Laser Ranging (SLR) technique ^Uj. Up to now, 
the only performed tests are due to Ciufolini and coworkers |llM12llT3lll4| . 
In such papers a confirmation of the existence of the Lense-Thirring effect as 
predicted by the Einstein's General Theory of Relativity is claimed with a 
total 20-25% (node-node-perigee LAGEOS-LAGEOS II linear combination) 
and 5—10% (node-node LAGEOS-LAGEOS II linear combination) accuracy, 
respectively, according to the adopted observable. 

In this paper we will analyze the latest results presented in |12[ I13[ I14j 
from a critical point of view in order to show that the total error budgets 
might have been underestimated. 

2 The Lense-Thirring effect on the orbit of a test 
particle and the strategy to measure it 

The gravitomagnetic field of a spinning mass of proper angular momentum 
J induces tiny secular precessions on the longitude of the ascending node O, 
and the argument of pericentre 2 uj of a test particle El EH El 

■ _ 2GJ . QGJcosi 

~ C 2 a 3(l-e 2 ) 3 / 2 ' ^ ~ " C 2 a 3(l_ e 2)3/2' W 

where G is the Newtonian constant of gravitation, c is the speed of light in 
vacuum, a, e and i are the semimajor axis, the eccentricity and the inclina- 
tion, respectively, of the test particle's orbit. See Figure ^ for the orbital 

1 According to K. Nordtvedt pQ, the multidecadal analysis of the Moon'orbit by means 
of the Lunar Laser Ranging (LLR) technique yields a comprehensive test of the various 
parts of order 0(c~ 2 ) of the post-Newtonian equation of motion. The existence of gravit- 
omagnetism as predicted by the Einstein's General Theory of Relativity would, then, be 
indirectly inferred from the high accuracy of the lunar orbital reconstruction. 

2 In their original paper Lense and Thirring use the longitude of pericentre zu — Q. + oj. 
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Figure 1: Orbital geometry for a motion around a central mass. Here L 
denotes the orbital angular momentum of the particle of mass m, J is the 
proper angular momentum of the central mass M, II denotes the pericentre 
position, / is the true anomaly of m, which is counted from II, f2, uj and i 
are the longitude of the ascending node, the argument of pericentre and the 
inclination of the orbit with respect to the inertial frame {x,y,z} and the 
azimuthal angle </> is the right ascension counted from the x axis. 

geometry of a Keplerian ellipse. In the terrestrial space environment the 
gravitomagnetic precessions are very small: for the spherically symmetric 
geodetic SLR LAGEOS satellites, whose orbital parameters are listed in Ta- 
ble Q they amount to a few tens of milliarcseconds per year (mas yr _1 in 
the following) 

The extraction of the Lense-Thirring precessions from the orbit data 
analysis is very difficult due to a host of competing classical orbital pertur- 
bations of gravitational [131131131201 and non-gravitational |2Tl l22l l23l \M 
1251 I26| origin which have various temporal signatures and are often quite 
larger than the relativistic signal of interest. The most insidious ones are the 

Table 1: Orbital parameters of the existing LAGEOS and LAGEOS II and 
of the proposed LARES and their Lense-Thirring node precessions. 



Satellite 


a (km) 


e 


i° 17 lt (mas yr 1 ) 


LAGEOS 


12270 


0.0045 


110 31 


LAGEOS II 


12163 


0.0135 


52.64 31.5 


LARES 


12270 


0.04 


70 31 
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perturbations which have the same temporal signature of the Lense-Thirring 
precessions 3 , i.e. secular trends. Indeed, whatever the length of the adopted 
observational time span T b s is, they cannot be fitted and removed from the 
time series without removing the relativistic signal as well. Then, it is of 
the utmost importance to assess as more accurately and reliably as possible 
their aliasing impact on the measurement of the Lense-Thirring effect. 

It turns out that the perigees of the LAGEOS-like satellites are severely 
affected by the non-gravitational perturbations, contrary to the nodes. More- 
over, since the non-conservative forces depend on the structure, the shape 
and the rotational status of the satellite their correct modelling is not a 
trivial task and, as we will see later, introduces large uncertainties in the 
correct assessment of the error budget in some of the performed gravitomag- 
netic tests. 

2.1 The gravitational error 

The even (I = 2,4,6...) zonal (m = 0) harmonic coefficients Ji of the mul- 
tipolar expansion of the Earth's gravitational potential, called geopotential, 
induce secular precessions 4 on the node, the perigee and the mean anomaly 
of any near-Earth artificial satellite |23| which, of course, depend only on 
its orbital configuration and are independent of its physical structure. Such 
aliasing effects are many orders of magnitude larger than the Lense-Thirring 
precessions; the precision with which the even zonal harmonics are known in 
the currently available Earth gravity models |2l ^ l2l ^Epi3Tl E ^ E ^ E ^ I3 ^ IHfi] 
would yield errors amounting to a significant fraction of the Lense-Thirring 
precessions or even larger. 

Even more dangerous are the perturbations induced by the secular vari- 
ations of the low degree even zonal harmonics Jg, I = 2, 4, 6 [371 I38| . Indeed, 
such perturbations grow quadratically in time if the shifts in mas are con- 
sidered and linearly in time if the rates in mas yr _1 are considered. Their 
impact on the orbital elements of the LAGEOS satellites have been worked 
out in 39 . It turns out that, by using the results of the errors induced 
by J2 would amount to 8%, 14% and 5.4% for the nodes of LAGEOS and 

3 Also the perturbations which grow quadratically in time are, of course, very danger- 
ous. Those induced by the secular variations of the even zonal harmonics of the Earth's 
geopotential fall in this category, as we will see in detail in Section r3.2.1l Harmonic time- 
dependent perturbations with periods longer than the observational time span may also 
be insidious because they would resemble superimposed linear trends 17 . 

4 Also the subtle non-gravitational Yarkovsky-Rubincam force, which is due to the 
interaction of the Earth's electromagnetic IR radiation with the physical structure of the 
LAGEOS satellites, induces secular effects on their nodes and perigees |23|. 
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LAGEOS II and the perigee of LAGEOS II, respectively, over an observa- 
tional time span T b s of just one year at 1— a level. This clearly shows that 
it would be impossible to analyze single orbital elements. 

The time-dependent harmonic perturbations HO] are less dan- 

gerous because if their periods are shorter than the adopted observational 
time span they can be fitted and removed from the time series. The most 
insidious tidal perturbation is that induced by the even zonal (055.565) con- 
stituent which has a period of 18.6 years and whose nominal impact on the 
orbital elements of the LAGEOS satellites amounts to thousands of mas 
[17| . However, it turns out that it does not affect the observables which 
have been adopted for the performed Lense-Thirring tests because its main 
component is of degree £ = 2 and order m = (see Section 12.1.2(1 . 

2.1.1 The butterfly configuration 

A possible way to cope with the even zonal gravitational perturbations is 
represented by the use of a pair of satellites in the so-called butterfly con- 
figuration. In it the orbital planes of the two satellites are shifted 180° 
apart, all other orbital parameters being equal. In this case the sum of the 
nodes [13111] fiW + C}(*+180°) and the difference of the perigees [32] 1331 133] 

£(<) _ ^(i+180°) 

would allow to cancel exactly out the bias due to all the 
even zonal harmonics of the geopotential 5 ; the gravitomagnetic precessions 
would, instead, add up. Indeed, the Lense-Thirring node precessions are 
independent of i, while the classical geopotential node precessions depend 
on a multiplicative factor cosi common to all degrees £ and on sums of 
even powers of sini which are different for the various degrees £ |18j : the 
Lense-Thirring perigee precessions depend on cos i while the classical geopo- 
tential perigee precessions depend on even powers of cosi and sini [T8] . 
The butterfly configuration cannot be realized with the present-day existing 
Earth artificial satellites. In [30] it was proposed to launch a LAGEOS-like 
satellite-the former LAGEOS III which later became LARES 6 [IHj-with the 
same orbital parameters of LAGEOS apart from the inclination which should 
be supplementary. That idea is still alive, although it has not yet been ap- 
proved by any national space agency or institution. Recently, it has been 
proposed to adopt the LARES orbital configuration for the OPTIS relativ- 
ity mission [461 1471 148j currently under examination by the German Space 

5 If totally passive satellites are considered, the difference of the perigees would yield 
less accurate results than the sum of the nodes |44|. 

6 The eccentricity of LARES would amount to glares = 0.04 in order to perform other 
tests of post-Newtonian gravity with the perigee. 
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Agency (DLR). It should be noted that with a LAGEOS-LARES/OPTIS 
butterfly configuration the difference of the perigees would not be a good 
observable because of the small eccentricity of the LAGEOS orbit. The 
exact cancellation of the classical precessions due to all the even zonal har- 
monics of geopotential could be achieved only if the LARES/OPTIS orbital 
parameters were exactly equal to those of LAGEOS. This would pose severe 
restrictions in term of the quality and, consequently, the cost of the rocket 
launcher to be used. Indeed, the realistically obtainable precision with the 
originally proposed sum of the nodes would be affected by the unavoidable 
departures of the LARES/OPTIS orbital parameters from their nominal 
values, at least to a certain extent. Also the fact that the eccentricities 
of LAGEOS and LARES/OPTIS would differ by one order of magnitude 
should be accounted for. Indeed, in these cases a residual aliasing effect 
due to all the even zonal harmonics J2, J4, J§... would still occur. These 
topics have been investigated in - A different observable involving also 
the nodes of LAGEOS and LAGEOS II has been proposed in [5U] for the 
LARES mission. It would dramatically reduce the dependence of the sys- 
tematic error due to the even zonal harmonics of the geopotential on the 
unavoidable orbital injection errors allowing also for a reduction of the costs 
of the mission. Moreover, as it will become clear later, also the impact of 
the Ji would be greatly reduced. On the contrary, the simple sum of the 
nodes would be affected by J2, J4, J§... which might not be negligible over a 
time span many years long. 

2.1.2 The linear combination approach 

The problem of reducing the impact of the mismodeling in the even zonal 
harmonics of the geopotential with the currently existing satellites can be 
coped in the following way. 

Let us suppose we have at our disposal N (N> 1) time series of the resid- 
uals of those Keplerian orbital elements which are affected by the geopoten- 
tial with secular precessions, i.e. the node and the perigee: let them be ^> A , 
A=LAGEOS, LAGEOS II, etc. Let us write explicitly down the expressions 
of the observed residuals of the rates of those elements Sip^ hs in terms of 
the Lense-Thirring effect ipjjp, of N-l mismodelled classical secular preces- 
sions ip^5Ji induced by those even zonal harmonics whose impact on the 
measurement of the gravitomagnetic effect is to be reduced and of the re- 
maining mismodelled phenomena A which affect the chosen orbital element 
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^dbs = ^lt/* + Yl ^ 6Ji + AA ' A = LAGEOS, LAGEOS II,... (2) 

N— 1 terms 

The parameter 7 fi is equal to 1 in the General Theory of Relativity and 
in Newtonian mechanics. The coefficients ip^ are defined as 

and have been explicitly worked out for the node and the perigee up to degree 
£ = 20 in jlUl I18j : they depend on some physical parameters of the central 
mass (GM and the mean equatorial radius R) and on the satellite's semi- 
major axis a, the eccentricity e and the inclination i. We can think about 
eq. ((21) as an algebraic nonhomogeneuous linear system of N equations in N 
unknowns which are /i and the N-l 5Jf. solving it with respect to \x allows to 
obtain a linear combination of orbital residuals which is independent of the 
chosen N-l even zonal harmonics. In general, the orbital elements employed 
are the nodes and the perigees and the even zonal harmonics cancelled are 
the first N-l low-degree ones. 

This approach is, in principle, very efficient in reducing the impact of the 
systematic error of gravitational origin because all the classical precessions 
induced by the static and time-dependent parts of the chosen N-l Ji do not 
affect the combination for the Lense-Thirring effect. Moreover, it is flexible 
because it can be applied to all satellites independently of their orbital con- 
figuration, contrary to the butterfly configuration in which the cancellation 
of the even zonal harmonics can be achieved only for supplementary orbital 
planes and identical orbital parameters. This approach has been adopted, 
among other things, in |5U1 147j in order to make the requirements on the 
LARES/OPTIS orbital configuration less stringent. Apart from the first 
orbital element which enters the combination with 1, the other elements are 
weighted by multiplicative coefficients Cj(a, e,i) ^ 1 which are built up with 
ip.e and, then, depend on the orbital elements of the considered satellites. 
Their magnitude is very important with respect to the non-gravitational per- 
turbations, which in general are not cancelled out by the outlined method, 
and to the other time-dependent perturbations of gravitational origin which 
are neither even nor zonal so that they affect the obtained combination as 
well. Values smaller than 1 for the q coefficients are, in general, prefer- 
able because they reduce the impact of such uncancelled perturbations. It 

7 It can be expressed in terms of the PPN 7 parameter |51| as /i = (1 + 7)/2. 
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is important to note that the order with which the orbital elements enter 
the combination is important: indeed, while the systematic error due to the 
even zonal harmonics of the geopotential remains unchanged if the orbital 
elements of a combination are exchanged, the coefficients q do change and, 
consequently, also the non-gravitational error. The best results are obtained 
by choosing the highest altitude satellite as first one and by inserting the 
other satellites in order of decreasing altitudes. 

This method was explicitly adopted for the first time in [52] with the 
nodes of the LAGEOS satellites and the perigee of LAGEOS II. The obtained 
combination is 

rrSLAGEOS , „ roLAGEOS II , „ r- LAGEOS II „ fin 9 ( A\ 

where c\ = 0.304, C2 = —0.350 and 60.2 is the slope, in mas yr^ 1 , of the 
expected gravitomagnetic linear trend. Eq. (|IJ) is insensitive to the first two 
even zonal harmonics J2 and J4. It has been used in when the level of 
accuracy of the JGM3 |2H] and EGM96 Earth gravity models, available 
at that time, imposed the cancellation of the first two even zonal harmonics, 
at least. 

In view of the great improvements in the Earth gravity field modelling 
with the CHAMP and, especially, GRACE |54j missions an extensive 
search for alternative combinations has been subsequently performed |55l 
l39| 156"! Ib7| 158] . In [221 EH! the following combination has been proposed 8 

^LAGEOS + ^LAGEOS II „ (g) 

where k\ = 0.546 and 48.2 is the slope, in mas yr , of the expected grav- 
itomagnetic linear trend. It has been adopted for the test performed in ^H] 
with the 2nd generation GRACE-only EIGEN-GRACE02S Earth gravity 
model |34j . Eq. © allows to cancel out the first even zonal harmonic J2. 

3 The performed Lense-Thirring tests with the 
LAGEOS satellites 

The only performed tests aimed at the detection of the Lense-Thirring pre- 
cessions of eq. (^Q) in the gravitational field of the Earth with the exist- 
ing LAGEOS satellites have been performed, up to now, by Ciufolini and 

8 The possibility of using only the nodes of the LAGEOS satellites in view of the im- 
provements in the Earth gravity models from GRACE has been propsed for the first time 
in |54|. although without quantitative details. 
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coworkers. They have used the node-node-perigee combination of eq. (@J) 
[111 112] and the node-node combination of eq. (J5J) ,13 . In these works it has 
often been claimed that "[...] the Lense-Thirring effect exists and its exper- 
imental value, [...], fully agrees with the prediction of general relativity. 11 

The main objections to the results presented in these works can be sum- 
marized as follows 

• The authors have not performed really robust and reliable tests e.g. 
by varying the length of the adopted observational time span, running 
backward and forward the initial epoch of the analysis, varying the 
secular rates of the even zonal harmonics in order to check their impact 
over different time spans, using different Earth gravity models in order 
to obtain a scatter plot of the obtained results. Instead, it seems that 
the authors, for a given data set, have always used from time to time 
those Earth gravity models which yielded just the closest results to 
what it is a priori expected from the General Theory of Relativity. 
Moreover, the impact of a priori 'imprint' effects of the Lense-Thirring 
signature itself on the coefficients of the Earth gravity models used 
may have driven the outcome of the performed tests just towards the 
expected result 

• The total error budget has probably been underestimated, especially 
the systematic error of gravitational origin. E.g., the impact of the sec- 
ular variations of the even zonal harmonics of the geopotential, which 
may become a very limiting factor over time spans many years long 
as those used, has not been addressed in a realistic, reliable and sat- 
isfactorily manner. Moreover, in some cases it has been incorrectly 
calculated by summing in quadrature the various contributions to the 
gravitational error (static even zonal harmonics, tides, secular vari- 
ations of the even zonal harmonics) which can, instead, hardly be 
considered as uncorrelated. Almost always 1— a results have been pre- 
sented without any explicit indication of this fact. 

• All the relevant works of other authors, in which many of these issues 
have been addressed, have always been consistently ignored 

• The node-node combination of eq. (JSJ) has been repeatedly and explic- 
itly presented as a proper own result of the authors with references to 
their works (ref. [6] of an d ref. 19 of ^3]) which, instead, have 
nothing to do with eq. © 
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3.1 The node-node-perigee tests 

The combination of eq. @ has been analyzed by using the EGM96 |29| 
Earth gravity model over 4 years in and over 7.3 years in [T2|. The 
claimed total error budget amounts to 20-25% over 4 years and to 20% over 
7.3 years. 

3.1.1 The gravitational error 

The impact of the remaining uncancelled even zonal harmonics of the geopo- 
tential Jq,Js, Jio, ... on eq. (jlj) has been estimated by Ciufolini and cowork- 
ers with the full covariance matrix of EGM96 in a root-sum-square calcula- 
tion. In and, six years later, in ^2] it is claimed to be < 13%. Apart 
from the fact that this is a 1 — a level estimate, in |59) . as later acknowledged 
in a number of papers [55J EEEI IM1 HSU 1551 158] . the use of the full covariance 
matrix of EGM96 has been questioned. Indeed, it has been noted that in 
the EGM96 solution the recovered even zonal harmonics are strongly recip- 
rocally correlated; it seems, e.g., that the 13% value for the systematic error 
due to geopotential is due to a lucky correlation between and J$ which 
are not cancelled by eq. (@J). The point is that, according to [5^j, nothing 
would assure that the covariance matrix of EGM96, which is based on a 
multi-year average that spans the 1970, 1980 and early 1990 decades, would 
reflect the true correlations between the even zonal harmonics during the 
particular time intervals of a few years adopted in the analyses by Ciufolini 
and coworkers. Then, a more conservative, although pessimistic, approach 
would be to consider the sum of the absolute values of the errors due to 
the single even zonal as representative of the systematic error induced by 
our uncertainty in the terrestrial gravitational field according to EGM96 
[23 EH]- In this case we would get a conservative upper bound of 83% (1-cr). 
If a root-sum-square calculation is performed by neglecting the correlations 
between the even zonals a 45% 1-cr error is obtained [Ts j 1551 139| 158]. 

3.1.2 The non-gravitational error 

Another important class of systematic errors is given by the non-gravitational 
perturbations which affect especially the perigee of LAGEOS II. The main 
problem is that it turned out that their interaction with the structure of 
LAGEOS II changes in time due to unpredictable modifications in the phys- 
ical properties of the LAGEOS II surface (orbital perturbations of radia- 
tive origin, e.g. the solar radiation pressure and the Earth albedo) and in 
the evolution of the spin dynamics of LAGEOS II (orbital perturbations of 
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thermal origin induced by the interaction of the electromagnetic radiation 
of solar and terrestrial origin with the physical structure of the satellites, in 
particular with their corner-cube retroreflectors) . Moreover, such tiny but 
insidious effects were not entirely modelled in the GEODYN II software at 
the time of the analysis of [111 I12| , so that it is not easy to correctly and 
reliably assess their impact on the total error budget of the measurement 
performed during that particular time span. According to the evaluations 
in 23] , the systematic error due to the non-gravitational perturbations over 
a time span of 7 years amounts to almost 28%. However, according to 
their impact on the measurement of the Lense-Thirring effect with the nodes 
of LAGEOS and LAGEOS II and the perigee of LAGEOS II is, in general, 
quite difficult to be reliably assessed. 

So, by adding quadratically the gravitational and non-gravitational er- 
rors of [23] we obtain for the systematic uncertainty <5^ s y stematlc 54% if we 
assume a 45% error due to geopotential. The sum of the absolute values of 
the errors due to gepotential added quadratically with the non-gravitational 
perturbations would yield a total systematic error of ^^ s y stematlc 88%. It 
must be noted that the latter estimate is rather similar to those released 
in [S3]. Note also that they are 1-cr evaluations. Moreover, it should be 
considered that the perigee of LAGEOS II is also sensitive to the eclipses 
effect on certain non-gravitational perturbations. Such features are, gener- 
ally, not accounted for in all such estimates. An attempt can be found in 
|21| in which the impact of the eclipses on the effect of the direct solar radi- 
ation pressure on the LAGEOS-LAGEOS II Lense-Thirring measurement 
has been evaluated: it should amount to almost 10% over an observational 
time span of 4 years. 

3.2 The node-node tests 

The situation might turn out to be even worse for the results presented 
for the node-node combination of eq. Such observable only cancels 

out the gravitational bias of the first even zonal harmonic J2, but has the 
great advantage of discarding the perigee of LAGEOS II and its insidious 
non-gravitational perturbations. 

It has explicitly been proposed for the first time in |56l I39j . although 
the possibility of using a LAGEOS-LAGEOS II node-only observable was 
presented for the first time in [SI] without quantitative details. In it 
seems that the author refers to it as a proper own result with the ref. [6], i.e. 
[TT] of the present paper. In ^3] the authors, instead, explicitly attribute it 
to themselves with ref. 19, i.e. [10] of the present paper. In [LT an unfair 
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use also of is present. 

3.2.1 The gravitational error 

In the node-node combination of eq. (JSJ) has been analyzed with the 
EIGEN2 [S21 and GGM01 [35] Earth gravity models over a time span of 
almost 10 years. 

In 39] the impact of the static part of the geopotential, according to the 
CHAMP-only EIGEN2 Earth gravity model, is evaluated as 18% (1-cr root- 
sum-square covariance calculation), 22% (1-cr root-sum-square calculation) 
and 37% (1-cr upper bound). Ciufolini only reports 18% obtained with 
the full covariance matrix of EIGEN2 for which the same remarks as for 
EGM96 holds. Moreover, he seems to ignore that EIGEN2 is only based on 
six months of data and that the released sigmas of the even zonal harmonics 
of low degree, which are the most relevant in this kind of analyses with 
the LAGEOS satellites, are rather optimistic, as pointed out in |32| and 
acknowledged in |39j . In regard to the GGM01 model, the covariance matrix 
was not publicly released. Cufolini correctly presents a 19% which is the 1—a 
upper bound obtained in [59*] . However, GGM01 is only based on 111 days 
of data. The author's claim "We conclude, using the Earth gravity model 
EIGEN-2S, that the Lense-Thirring effect exists and its experimental value, 
fi = 0.98 ± 0.18, fully agrees with the prediction of general relativity" seems 
too optimistic. 

In 13: the authors use the 2nd generation GRACE-only EIGEN-GRACE02S 
Earth gravity model |34j . Also in this case the full covariance matrix was 
not available and the authors correctly report a systematic error due to the 
even zonal harmonics of 3% (root-sum-square calculation) and 4% (upper 
bound) at 1-a level. The same results have been obtained in |60j . 

Major problems may arise when the authors show their a priori error 
analysis for the time-dependent gravitational perturbations (solar and lunar 
Earth tides, secular trends in the even zonal harmonics of the Earth's field 
and other periodic variations in the Earth's harmonics). Indeed, in |13j 
they claim that, over an observational time span of 11 years, their impact 
would be 2%. This evaluation is based on ref. 30 of ^3] which refers to the 
WEBER-SAT/LARES INFN study; it has nothing to do with the present 
node-only LAGEOS-LAGEOS II combination. 
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3.2.2 The impact of the secular rates of the uncancelled even 
zonal harmonics 

Moreover, this estimate may turn out to be optimistic because of the secular 
variations of the even zonal harmonics 9 Jg. Indeed, eq. © allows to cancel 
out J2, but, in principle, is sensitive to J4, Jq,-.-, as pointed out in [SO] . 
The uncertainties in the Ji are still quite large: e.g., according to Table 1 of 
|38| . there is not yet even full consensus on the sign of such rates. On the 
other hand, their impact on the Lense-Thirring measurement grow linearly 
in time 10 . Indeed, the mismodelled shift, in mas, of eq. (JHJ due to the 
secular variations of the uncancelled even zonal harmonics can be written 
as 

£ (^ AGEOS + kl n¥ GEOS n ) 6 ~yTl s , (6) 

1=2 

where the coefficients Clj are <9f2 c i ass /dJ£ and have explicitly been calculated 
up to degree £ = 20 in ^IJ^]. It must be divided by the gravitomagnetic 
shift, in mas, of eq. ® over the same observational time span 

(O^GEOS + ^LAGEOS ifj ^ = ^ mag yf -l ^ (?) 

By assuming 5J 4 = 0.6 x 10~ n yr" 1 and 6J G = 0.5 x 10~ n yr" 1 [HHJ , it 
turns out that the percent error on the combination eq. © grows linearly 
with T Q bs and would amount to 1% over one year at 1 — a level. This means 
that, over 11 years, their impact might range from 11% (1-cr) to 33% (3- 
a). These evaluations hold on the assumption that J4 and Jq retain their 
constant sign over the adopted observational time span. At present, there 
is, in fact, evidence for an inversion only in J2 since 1998 |38j . but eq. © 
is insensitive to it. In Jl] some explanations about the procedure followed 
in the analysis of |13j can be found. Basically, they follows 

• The time series of the combined residuals of the nodes of LAGEOS and 
LAGEOS II, built up without the even zonal rates in the background 
reference model for the Earth gravity field, has been fitted with some 
harmonic signals, a linear trend and a parabolic signal from which a 

9 The problem of the secular variations of the even zonal harmonics in post-Newtonian 
tests of gravity with LAGEOS satellites has been quantitatively addressed for the first 
time in |61|. In regard to the Lense-Thirring measurement with eq. ©, it has been, 
perhaps, misunderstood in 39 . 

1 For a possible alternative combination which would cancel out the first three even 
zonal harmonics along with their temporal variations see |57l I58| . 
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jf ~ 1.5 xlO -11 yr" 1 has been determined. Then, the so obtained jf 
would have been inserted in the dynamical force models of the orbital 
processor in order to re-analyze the same data yielding a certain value 
of jj, which we conventionally define as fi ™a • 

• The time series of the combined residuals of the nodes of LAGEOS 
and LAGEOS II, built up with the default value of the even zonal 
rate jf = 1.41 x KT 11 yr" 1 in the EIGEN-GRACE02S model, has 
been fitted, again, with some harmonic signals, a linear trend and a 
parabolic signal yielding a certain value of \i which we conventionally 

define aS fl jdefault • 

• The so obtained values \i j eR and ^default have been compared by find- 
ing a 1% variation only. The author of writes " [...] using the value 
jEffective ^ Lg # 10 _n for thg LAGEOS satellites that we obtained 
from fitting the combined residuals [...] resulted in a change of the 
measured value of frame-dragging by about 1% only with respect to 
the case of using j 4 = 1.41 * 10~ n [...]" 

3.2.3 Some possible criticisms to the outlined approach 

The value jf ~ 1.5 x 1CP 11 yr _1 measured with the combination of eq. ® 
is affected not only by Jq and the other higher degree even zonal harmonics 
but, more importantly, by the Lense-Thirring signature itself. Indeed, the 
combination of eq. © is designed in order to only disentangle J2 and the 
Lense-Thirring effect. So, it is not admissible to use the so obtained Jf , 
which is coupled by construction to the Lense-Thirring effect, in order to 
reliably and correctly measure the Lense-Thirring effect itself. A combina- 
tion suitably designed in order to measure J4 independently of J2, Jq and 
the General Theory of Relativity has been proposed in |H2] 

3.2.4 Some suggestions to improve the reliability and the robust- 
ness of the discussed test 

• The values of Ji to be used in the dynamical force models of the or- 
bital processors should have been determined in a Lense-Thirring-free 
fashion in order to avoid 'imprinting' effects on the outcome of the 
analysis. Indeed, they could easily drive it just towards the expected 
result. In regard to this point, suitably designed linear combinations 
of the residuals of some orbital elements of the existing SLR satel- 
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lites could be used in order to disentangle, e.g., J2, J4, Jq and the 
gravitomagnetic force 

• For a given observational time span and a given background refer- 
ence Earth gravity model, two time series of the combined residuals, 
built up with and without the aforementioned Ji GTR-free values in 
the background reference models, should be analyzed The difference 
between the so obtained /j, parameters can, then, be evaluated 

• The aforementioned procedure should be repeated for various obser- 
vational time spans in order to have a more robust assessment of the 
importance of the secular variations of the even zonal harmonics in the 
proposed measurement. Indeed, their systematic error on the Lense- 
Thirring effect grows linearly 

• The aforementioned procedures should be repeated for various Earth 
gravity models in order to have a scatter plot 

3.2.5 On the correct evaluation of the systematic error of gravi- 
tational origin 

Another controversial point is that it is unlikely that the various errors 
of gravitational origin can be summed in a root-sum-square way because of 
the unavoidable correlations between the various phenomena of gravitational 
origin. Indeed, in the adopted EIGEN-GRACE02S model J4 and Jq are not 
solved for: J2 and J4 have been assumed fixed, while Jq is absent. This 
means that the recovered are affected by such rates. It would, then, be 
more conservative to linearly add the errors induced by them. In this case, 
the ( jf ] - J £ ) error would range from 15% (4%+ll%) at 1 — a level to 45% 
(12%+33%) at 3 — (j level over a 11-years long observational time span 11 . 
The so obtained global gravitational error can be added in quadrature to 
the non-gravitational error. Even by assuming the 2% authors' estimate of 
the time-dependent part of the gravitational error, the upper bound errors 
would be y/(4 + 2) 2 + 2 2 % = 6% at 1-ct, ^/(8 + 4) 2 + 4 2 % = 13% at 2-a 
and \J (12 + 6) 2 + 6 2 % = 19% at 3-cr. Instead, in [IB] the authors add in 
quadrature the doubled error due to the static part of the geopotential (the 
2x4% value obtained from the sum of the individual error terms), their 
perhaps optimistic evaluation of the error due to the time dependent part of 
the geopotential and the non-gravitational error getting \/8 2 + 4 2 + 4 2 % = 
10% at 2-<T. On the other hand, when they calculate the 3-a upper bound it 

11 The evaluations of eq. J|jJ have been used. 
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seems that they triple the 3% error due to the static part of the geopotential 
obtained with a root-sum-square calculation and add it in quadrature to the 
other (not tripled) errors getting V9 2 + 2 2 + 2 2 % < 10% at 3-<r. 

3.2.6 The a priori 'memory' effect of the Lense-Thirring signa- 
ture on the adopted Earth gravity model 

The following point seems also worthy of discussion. The recovered values 
of the even zonal harmonics in the GRACE models like EIGEN-GRACE02S 
retain, in principle, a sort of 'imprint' or 'memory' of the General Theory of 
Relativity, which, in fact, has not been modelled in the currently released 
GFZ-GRACE-based models (F. Flechtner, GFZ team, private communica- 
tion, 2004). This is certainly true for the Earth gravity models of the pre- 
CHAMP and GRACE era obtained from multidecadal laser ranging to the 
geodetic satellites of LAGEOS type, as discussed in [HS! EI2 EE] • This fea- 
ture is also valid for, e.g. GRACE which recovers the low degree even zonal 
harmonics from the tracking of both satellites by GPS and the medium-high 
degree geopotential coefficients from the observed intersatellite distance vari- 
ations. Indeed, from 6i\ it can be noted that the variation equations for the 
Satellite-to-Satellite Tracking (SST) range Ap and range rate Ap of GRACE 
can be written in terms of the in-plane radial and, especially, along-track 
components R, T, Vr,Vt of the position and velocity vectors, respectively. 
In turns, they can be expressed as functions of the perturbations in all the 
six Keplerian orbital elements (see (10)-(11), (A4)-(A6), (A14)-(A16) and 
(A28)-(A30) of |64j ). Now, the gravitomagnetic off-diagonal components 
of the spacetime metric also induce short-periodic 1 cycle per revolution (1 
cpr) effects |Ul 1631 Ho] on all the Keplerian orbital elements, apart from the 
secular trends on the node and the pericentre. This means that there is also 
a Lense-Thirring signature in all the other typical satellite and intersatellite 
observables like ranges and range-rates. It is likely that it mainly affects 
the low-degree even zonal harmonics. In this case, it might happen that the 
results of H3U3], obtained with the J 2 -free node-node LAGEOS-LAGEOS 
II combination, have been driven close to the expected result, at least to a 
certain extent, just by this a priori Lense-Thirring 'imprint' on the Earth 
gravity model used. 

3.2.7 The LARES mission 

Finally, it is hard to understand why the authors of [Hj very often refer to the 
LAGEOS-LARES proposed experiment and to the related simulations and 
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error budgets. It is rather confusing and misleading. As already explained, 
the LAGEOS-LAGEOS II combination of eq. (JSJ) is, by construction, de- 
signed in order to exactly cancel out the J2 term with an approach which can 
be applied to any orbital configuration given a pair of satellites in different 
orbits or a pair of different Keplerian orbital elements of the same satellite 
(see Section [2.1.2(1 . On the contrary, the observable originally proposed for 
the LAGEOS-LARES mission is the simple sum of their nodes. As pointed 
out in Section [2.1.11 if the orbital parameters of LARES, quoted in Table 
n were exactly equal to their nominal values, all the even zonal harmonics 
would be exactly cancelled out. Instead, the sum of the nodes would be 
affected, to a certain extent, by the whole range of the even zonal harmonics 
of the geopotential due to unavoidable departures from the LARES nominal 
configuration because of the orbital injection errors and mission design (the 
eccentricity of LARES would be one order of magnitude larger than that 
of LAGEOS), i.e. the coefficients of the classical nodal precessions would 
not be exactly equal and opposite f^AGEOS / -fjLARES for t = 2, 4, 6, 8, .... 
In ^3] the combination of the nodes of LAGEOS and LAGEOS II of eq. 

is presented as if it is only slightly different with respect to the sum of 
the nodes of the originally proposed LAGEOS-LARES configuration, apart 
from a 18° offset in the inclination of LAGEOS II with respect to LARES. 
The differences in the eccentricities and the semimajor axes, which do play 
a role have been neglected. 

4 Conclusions 

In this paper we have performed a detailed critical analysis of the reliability 
and robustness of the so far performed tests aimed at the detection of the 
Lense-Thirring effect in the gravitational field of the Earth with the existing 
or proposed LAGEOS satellites. 

We can summarize our conclusions as follows 

• In regard to the node-node-perigee LAGEOS-LAGEOS II combina- 
tion, still presented in recent works as J2j, the claimed 20 — 25% total 
accuracy obtained with the EGM96 Earth gravity model seems to be 
too optimistic mainly because of an underestimation of the systematic 
error due to the geopotential. 

• In regard to the node-node LAGEOS-LAGEOS II combination of eq. 
©, extensive and thorough tests are required in order to check the 
impact of the secular variations of the even zonal harmonics of the 
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geopotential which may represent a very limiting factor over time spans 
many years long. Such tests should be conducted by varying the data 
sets, shifting backward and forward the initial epoch of the analyses, 
varying the magnitudes of the Ji in the force models for a given ob- 
servational time span and the length of the time span for given values 
of Ji. If it will turn out that J4, Jq, ... do affect eq. only improve- 
ments in our knowledge of the Jg will allow to use this observable in a 
truly reliable and confident way; a better knowledge only of the static 
part of the geopotential would not be sufficient. Moreover, also the 
problem of the Lense-Thirring 'memory' in the adopted gravity model 
should be accounted for. 

• Alternative combinations involving the use of existing SLR targets 
other than the LAGEOS satellites should be analyzed. The most 
promising combination involves the nodes of LAGEOS, LAGEOS II, 
Ajisai and Jason-1 [23 EH]- It cancels out the first three even zonal 
harmonics J2, J4, Jq, along with their temporal variations, at the price 
of introducing the relatively huge non-gravitational perturbations on 
Jason-1 which, however, should have a time-dependent harmonic sig- 
nature with short periodicities. Moreover, in this case the possible im- 
pact of the Lense-Thirring 'imprint' would be greatly reduced. Accord- 
ing to the recently released combined CHAMP+GRACE+terrestrial 
gravimetry/altimetry EIGEN-CG01C Earth gravity model (3^], the 
systematic error due to the remaining even zonal harmonics would 
amount to 0.7 % (root-sum-square calculation) and 1.6% (upper bound) 
at 1— a |57| . The possibility of getting long time series of the Jason's 
node should be seriously investigated with real data tests 

• The launch of a third LAGEOS-like satellite would allow for a really 
robust and confident measurement of the Lense-Thirring effect. In- 
deed, it would be possible to combine its node-and, perhaps, also its 
perigee if it will be built up in such a way to sufficiently reduce the 
impact of the non-gravitational perturbations-with the nodes of the 
other currently existing LAGEOS satellites in order to cancel out J2, 
J4 and, perhaps, Jq along with their temporal variations. The great 
improvements in our knowledge of the terrestrial gravitational field 
would allow to abandon the very tight and stringent requirements of 
the originally proposed butterfly configuration. However, since the 
LARES data should be combined with those of the existing LAGEOS 
satellites, the limit of the total available accuracy is mainly set by the 
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impact of the non-gravitational perturbations, i.e. it is of the order 
of 1%. The systematic error of gravitational origin could be kept well 
below this value. 
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